arXiv:1504.05700vl [astro-ph.SR] 22 Apr 2015 


MAGNETOHYDRODYNAMIC SHOCKS IN AND ABOVE 
POST-FLARE LOOPS: TWO-DIMENSIONAL SIMULATION 

AND A SIMPLIFIED MODEL 


Shinsuke Takasao 1 , Takuma Matsumoto 2 , Naoki Nakamura 1 , and Kazunari Shibata 


1 


takasao@kwasan.kyoto-u.ac.jp 


Received 


accepted 


1 Kwasan and Hida Observatories, Kyoto University, Yamashina, Kyoto 607-8471, Japan 

institute of Space and Astronautical Science (ISAS), Japan Aerospace Exploration 
Agency (JAXA), 3-1-1 Yoshinodai, Sagamihara, Kanagawa 252-5210, Japan 





2 


ABSTRACT 


Solar flares are an explosive phenomenon, where super-sonic flows and shocks 
are expected in and above the post-flare loops. To understand the dynamics 
of post-flare loops, a two-dimensional magnetohydrodynamic (2D MHD) simu¬ 
lation of a solar flare has been carried out. We found new shock structures in 
and ab ove t he p ost- flare loops, which were not resolved in the previous work by 


Yoko vama fe Shi bata fl2001h . To study the dynamics of flows along the recon¬ 


nected magnetic field, kinematics and energetics of the plasma are investigated 
along selected field lines. It is found that shocks are crucial to determine the 
thermal and flow structures in the post-flare loops. On the basis of the 2D MHD 
simulation, we have developed a new post-flare loop model which we call the 
pseudo-2D MHD model. The model is based on the ID MHD equations, where 
all the variables depend on one space dimension and all the three components of 
the magnetic and velocity fields are considered. Our pseudo-2D model includes 
many features of the multi-dimensional MHD processes related to magnetic re¬ 
connection (particularly MHD shocks), which the previous ID hydrodynamic 
models are not able to include. We compare the shock formation and energetics 
of a specific field line in the 2D calculation with those in our pseudo-2D MHD 
model, and we found that they give similar results. This model will allow us to 
study the evolution of the post-flare loops in a wide parameter space without ex¬ 
pensive computational cost and without neglecting important physics associated 
with magnetic reconnection. 


Subject headings: Sun: corona — magnetic fields — stars: flare — magnetic 
reconnection — Sun: flares — Sun: oscillations — shock waves — waves - 
magnetohydrodynamics (MHD) 
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Introduction 


Solar flares are a phenomenon in which a huge amount of magnetic energy stored in the 
coronal held is rapidly released by magnetic reconnection. Where magnetic reconnection 
is a physical process in which a magnetic held in a highly conducting plasma changes its 
connectivity due to finite resistivity. The rapid energy release is accompanied by various 
plasma dynamics, like shock waves and the chromospheri c ev a porat ion. For more details 


about solar hares, the reader is referred to 


Shi bata & Ma gara (120111) . 


Flares similar to solar_ hares have 
astrophysical objects (e.g. Koyama et ah 


jeen observed from many stars and other 


1994; 

Tsuboi et al. 

1998; 

Giidel 

2004) 


of those hares and solar hares have many common features, solar hare models have been 


applied to those astrop 


Machida fe Matsnmoto 


rysica. 


,2003 


hares (e.g. 


H 


Masada et al. 


avashi et ah 


1996 


Shibata & Yokovama 


2002 


2010). Note that there should be differences in 


the plasma parameters between those hares and solar hares, like the plasma density, Alfven 
speed, and system size. Since these values determine the time scales of different processes, 
this can make the evolution of different hares distinct. Therefore, to comprehensively 
understand the hare physics in the universe, we need to explore various hares in a wide 
range of the plasma parameters. 

There are many attempts to model solar and stellar hare loops. One-dimensional (ID) 
hydrodynamic models have been widely developed to study the thermal evolution and hows 
in a single loop. Depending on the assumed main energy transfer p rocess, t he models can be 


categ orized into two branches: conduction-heating models 


1993) and clectron-beam-heating models (e.g. 


Fisher et ah 


e.g. 


1985 


agai 


1980 


Peres fe Reale 


Mariska et ah 


general, the results of both models agree with observed emissions well. 


Hori et ah 


1989). In 


(11997 


1998 ) developed a pseudo-2D loop model based on ID hydrod ynamic calculations, and 


gave a simple description of observed soft X-ray emissions (see 


Warren (120061 ) for further 
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development of this kind of models). 


The first 2D magnetohydrodynamic (MHD) simulation of a solar fl are based on a 


reconnection model with the heat conduction was performed by 


Yo kovama fe Shi bata 


(120011 ). and they studied the thermal evolution in 


that determines the flare temperature. 


Shiota et ah 


he post-flare loops and the physics 


reconnection with the heat conduction to study the coronal 
giant arcade, where they do not include the chromosphere. 


(120051 ) performed MHD simulations of 
mass ejections and associatec 


Miyagoshi & Yokovama (2001) 


carried out a MHD simulation of chromospheric evaporation jets formed as a result of 


Longcope et al. 

(2009 

) and 

Longcope 

(2014) 


dynamics of a reconnected field line is considered. 

Shocks in and above the post-flare loops can be important for both the non-thermal 
particle acceleration and evolution of the thermal structure. The above-the-looptop hard X- 
ray source found by 


Masuda et al. 


mechanism (for recent progress, see 


(11995) has posec 


Krucker et al. 


prob lems abou t the electron acce leration 


(e.g. 


2010) and 


Oka et al. 


(20151)). One 


of the promising scenarios is the acceleratio n by fast shocks which are expected to 


De 


formed above the soft X- ray post-flare loops dSomov fe Kosusi 


Tanuma fe Shibata 


19971 : 


Tsuneta fe Naito 


20051 ) . Recently acceleration from contracti ng plasmoids (Drake et al 


2006) and acceleration in plasmoids interacting with fast shocks (Nishizuka & Shibata 


1998 


2013) 


have been also proposed. It has been argued that 

a high-density 

region can be formed as a 

result of the shock-shock interaction at the top 

Lori et al. 

1997 

) and by the compression 

at termination shocks ( 

Yokovama & Shibata 

2001 

L). 

Zenitani & Mivoshi 

(2011 

) carefully 


analyze shocks in a plasmoid in a ideal 2D MHD simulation. Despite the importance, the 
shock formation in and above the post-flare loops has not been yet investigated in detail. 


How post-flare loops evolve in many astrophysical systems is one of our central 
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interests. Since reconnection produces hot and high-speed plasma flows, it is expected 
that the thermal structure should be determined by a coupling among the plasma flows, 
shocks, and heat conduction. The heat conduction is essential to drive the hot high-speed 
evaporation flows. However, the calculation of the plasma dynamics with the heat 
conduction effects is computationally expensive. This makes the extensive parameter 
survey in the multi-dimensional MHD simulations difficult. The development of a simplified 
post-flare loop model based on a reconnection model is therefore desired. 


In this paper, we study the shock formation and evolution of the thermal structure in 
and above the post-flare loops using MHD simulations. To understand the shock formation 
in 2D systems, a 2D MHD simulation of a solar flare has been carried out. We found 
new shock structures in and above t h e pos t-flare loops, which were not well resolved in 


the previous work by 


Yokovama & Shibata ( 20011 ). To study the dynamics of flows along 


the reconnected magnetic field in detail, kinematics and energetics of the plasma are 
investigated along selected field lines. It is found that shocks in the post-flare loops are 
crucial to determine the thermal and flow structures in the post-flare loops. On the basis 
of the results of the 2D MHD simulation, we have developed a new post-flare loop model 
which we call the pseudo-2D MHD model. The model is based on the ID MHD equations, 
and has been developed to model a 3D reconnected field line. Through a comparison, it is 
found that the shock formation and thermal evolution in the pseudo-2D MHD model are 
similar to those in the 2D MHD model. 


This paper is structured as follows. In Section [21 the 2D MHD model of a solar flare 
is introduced. In Section [3l numerical results of the 2D MHD simulation are detailed. 

In Section [H on the basis of the 2D MHD simulation, we develop our pseudo-2D MHD 
model of a post-flare loop. In Section [5l numerical results of the pseudo-2D MHD model 
are introduced and are compared with the 2D MHD simulation. Section [6] contains our 
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summary and discussion. 


2D MHD Model of a Solar Flare 


2.1. Assumptions and Basic Equations 


We performed a 2D MHD simulation of a solar flare similar to 


Yokovama fc Shibata 


u200ll) simulations. We take a rectangular calculation domain in the x-y plane. The MHD 
equations in the following form are solved: 


| + v. M =o, 


dpv 

~df 


pvv+ [p + 


\B\‘ 
8ir 

dB 


BB 


= 0 , 


dt 


de 


v e + p + 


w 

87T 


-j—B(v ■ B) + x B + F c 

471 ' 47T 


47T 

cV x E = 0, 
= 0, 


F c = -k 0 T 5 / 2 V||T, 

kB rji 

P = — pT, 
m 

E = - v x Bj 

P 

e =-- + -pv + 


7 


8tt’ 

J = —V x B, 

4:7T 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 

(9) 


where v = {v x ,v y ,v z ) and B = (B x ,B y ,B z ). k 0 ( 10 6 in cgs units) is the coefficient of 

the Spitzer-type conductivity, 7 is the specific heat ratio (5/3 is used in this study), m 
is the mean particle mass and fc# is the Boltzmann constant. // is the electric resistivity. 
F c is the conduction flux. The normalization units of our simulations are summarized in 
Table [Q The radiative cooling time is expected to be longer than the dynamical time in 
the post-flare loops. Since we mainly focus on the dynamical processes such as the shock 
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formation and flows, we neglect the radiative cooling. 


2.2. Initial and Boundary Conditions 


Our model is identical to the 


Yokovama & Shibata (120011 ) model, except for the domain 


size and boundary conditions. The initial condition is same as their model. The initial and 
boundary conditions are as follows. The domain is 0 < x < x max pD and 0 < y < y ma x,2D , 
where x maX) 2D= 10 and y m ax,2D = 20. A schematic diagram of the initial and boundary 
conditions are shown in Figure [U The initial magnetic field is assumed to be a force-free 
field and is given by 


B x (x, y) = 0, 

(10) 

B y (x,y) = Bq tanh (x/w), 

(11) 

B z (x, y) = B 0 / cosh (x/w), 

(12) 


where w = 0.5 is the width of the initial current sheet. The gas pressure is assumed to be 
uniform ( p 0 ). The density distribution is described as 

p(x, y) = p chr + (Pcor ~ Pchr (tanh [{y - Iitr)/w T r,2d\ + 1), (13) 

where p cor and p c h r are respectively the densities in the corona and chromosphere, and 
hxR and wtr, 2 D are respectively the height and the width of the transition region between 
the corona and chromosphere. h TR and Wtr^d are set to 1 and 0.2, respectively. The 
chromosphere is modeled as a dense and cool plasma. For simplicity, in this paper p c h r is 
set to 10 5 p cor . T cor = (m/kB)(po/Pcor) is the initial coronal temperature. 

To allow the magnetic field to reconnect, we impose a localized resistivity in the form 
of 


y(x,y) = po exp [-(r/w v ) 2 ], 


(14) 





where r = \Jx 2 + (y — and w v — 1. We localize the resistivity and £x it wi 


realize a fast and quasi-steady magnetic reconnection with a single X-point (jUgai 


1 time to 


1992). 


All the boundaries are symmetric. At the boundary at x = 0 the sign of B v is reversed. 


The numerical sc 


by 


Mivoshi & Kusano 


lerne is based on a Harten-Lax-van Leer (HLL) scheme developed 
(120051) . HLLD (”D” stands for Discontinuities), which is a shock¬ 


capturing scheme. It has the second-order accuracy in space and time. The heat conduction 


term is solved with an implicit method similar to 


Yokovama fc Shibata (20011) method to 


reduce the calculation time. We modify their method to more accurately calculate the heat 
conduction flux. A detailed description of the implicit method is given in Appendix |A] The 
calculated domain is resolved with 800 x 1000 grids. 


3. Numerical Results of 2D MHD Model 

3.1. Overview of Evolution 

Figure [2] displays an overview of the time evolution of the 2D MHD simulation. The 
region where x < 0 is also shown only for visual inspection. Reconnection takes place due 
to the localized resistivity and the Alfvenic collimated reconnection outflow is produced. 
The reconnected field lines pile up in the outflow region to form the growing loop system 
at the base of the corona. The temperature structure is smooth along the magnetic held 
owing to the heat conduction. The heat in the hot outflow is carried to the chromosphere 
along the magnetic held due to the heat conduction. As a result, the upper chromosphere is 
heated up and rapidly expands, leading to the formation of the hot dense uphows called the 
chromospheric evaporation. The evaporated plasma finally fills the reconnected held lines 
and forms the hot dense post-hare loops. We note that the weak disturbance which starts 
to propagate at the beginning of the simulation is generated because the initial condition is 
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not in the thermal equilibrium between the chromosphere and corona. We confirmed that 
its effect is negligible for the dynamics of the post-flare loops. 


One-dimensional plots parallel to the x —axis across y = 10 (i.e. across the reconnection 
outflow) are shown in Figure [3l It is shown that a Petschek-type reconnection is established 
because of the localized resistivity, where two slow shocks emanate from the reconnection 
region. The slow shocks can be discerned as a pair of the discontinuities in Figure [3j 


We note that due 


(Forbes et al. 


1989 


;o the heat conduction the s' 


Yokovama & Shibata 


19971; 


ow shocks 


Decome isothermal slow shocks 


Chen et al. 


1999; 


Seaton & Forbes 


2009). 


The thermal structure of the post-flare loops is determined by a complex coupling 
among the plasma flows, shocks, and heat conduction. When the reconnection outflow 
collides with the loop system below, two oblique fast-mode shocks and sometimes a 
horizontal fast-mode shock (Mach disk) are formed at the top (see Figure [4]), which 
indicates that the termination shock is a combination of two or three fast-mode shocks. 
Note that most of the regions where V • v/C s takes large negative values are slow or 
fast shocks. We found that the strength of the termination shock shows a quasi-periodic 
oscillation. We also found the shock reflection and Mach reflection in a concave magnetic 
field at the top (Figures [4] and |6|). The relationship between the oscillation and the flow 
structure at the top will be discussed in our future papers. 


The high-pressure plasma at the loop-top expands along the magnetic field to the 
foot-points, generating the high-speed downflows. Then the evaporation flows collide with 


the downf 
named by 


ows, forming the den se regions in the post-flare loops (” humps” in Figure El 


Yokovama fe Shibatal (120011) ). After the collision, the fronts of the evaporation 


flows and downflows becomes slow shocks (see V • v/C s maps). The pair of the upward 
slow shocks finally collide with each other at th e apex, forming another dense region (see 


Figure [6]). Note that the ’’blob” named by 


Yokovama & Shibata (2001J) is different from 
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this high-density region. The blob is the high-density region in a concave magnetic held at 
the top. 


Figure [7] displays the slow shocks in the post-hare loops mentioned above. The 
entropy is discontinuous at the two pairs of the discontinuities in the density map, but 
the temperature is smooth along the magnetic held, indicating that they are isothermal 
slow shocks. The plasma ft in a large fraction of the post-hare loops, as well as in the 
outhow region, is larger than unity (the contour indicates the level where ft = 1), meaning 
that the low-/3 approximation cannot be applied to the reconnected magnetic held. This is 
a consequence of the shock heating and compression. The high-pressure post-hare loops 
are co nfined by the surrounding law-ft plasma. This is consistent with the observation by 


McKenzie (120131 ): they concluded that plasma ft is of the order of unity or larger than unity 


in the supra-arcade plasma in two hares analyzed. 


3.2. Dynamics and Energetics along a Specific Field Line 

We have seen the two-dimensional evolution of a solar hare. To study the dynamics 
of hows along the reconnected magnetic held in detail, kinematics and energetics of the 
plasma are investigated along selected held lines. We performed the same analysis for other 
held lines and confirmed that they give similar results. 

We picked up a magnetic held line and measured the physical quantities along it. The 
held line used in the analysis is shown in Figure [8] Figure |9] displays the time-distance 
diagrams of the density, temperature, un = \(v ■ B)B/B 2 \ and V • v/C s along the held line 
whose foot-point is located at (x,y) = (1.2,0), where C s is the local sound speed. The 
distance is measured along the held line. Before the held line reconnects, the starting point 
of the distance is the intersection point between the held line and the top boundary. After 
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reconnection, the starting point is the apex of the closed loop. The sign of vii is defined as 
positive and negative when a plasma flow travels to the apex and foot-points, respectively. 

In Figure [9] we can discern the field line shrinking after reconnection. The heat released 
at the Petschek slow shocks is transferred to the chromosphere (see the temperature in 
Figure E]), leading to the generation of the chromospheric evaporation. The evaporation 
flows are seen as upflows from the chromosphere (see p and ny maps). The reconnection 
outflow is significantly decelerated at £ ~ 12. This is seen as the enhancement of the density 
and temperature and negative V • v (compression), indicating that the kinetic energy of the 
outflow is converted to the thermal energy. After the termination, the high-pressure plasma 
at the top expands along the magnetic held, forming the downflows. The downflows collide 
with the chromospheric evaporation hows, leading to the formation of the humps (see also 
Figure [5]). After the collision, the fronts of the evaporation hows and downflows become 
steep and finally become slow shocks. The pair of the upcoming shocks finally cross each 
other at the apex at £ ~ 16, forming another high-density region (see also Figure [6]). 


We found slow shocks propagating along the held line, but there is no prominent shocks 


nor waves propagating b ack a nc 


conduction (e.g. 


Ofman & Wang 


forth from end to end. They are damped by the heat 


200211 . Also note that the propagation speed of the shocks 


is largely decelerated by the evaporation how (i.e. Doppler effect), which significantly 
affects the traveling time of the slow-mode waves/shocks (the local sound speed is ~ 1.5, 
but the propagation speed is ~ 0.6). That is, a simple estimation by travel,slow ~ L/C s is 
not a good approximation of the traveling time, where L and C s is the loop length and the 
sound speed in the loop, respectively. We observe no prominent signature of the standing 


slow-mode waves in the ca 
see e.g. 


Nakariakov et a! 


culated time range (as for acoustic waves in the post-hare loops, 


(2004)). 


Figure [10] shows the time evolution of the total, magnetic, internal (thermal) and 
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kinetic energies integrated along two specific field lines. The top and bottom panels are for 
the field lines which are rooted at (x,y) = (1.2,0) and (1.8,0), respectively. Considering 
that the cross-sectional area of the flux tube is inversely proportional to the magnetic field 
strength, we integrate the energies along a field line as follows: 


Emag 


(15) 

Eint 


(16) 

Ekin 


(17) 

Etot = 

Emag Ei n t ~h Ekin' 

(18) 


In this paper, all the energies are normalized by the initial total energy. After the field 
lines reconnect (t ~ 8), the magnetic energy is rapidly converted to the internal and kinetic 
energies. As for the field line rooted at (x,y) = (1.2,0), the maximum value of the kinetic 
energy is 0.35. Let us define the time when E kin becomes E kin rnax as t peak . If we compare the 
total, magnetic and internal energies at t — 0 and those at t — t peak + 10, dE tot = E tot (t = 
tpeak + 10) — E to t(t = 0) ~ 0.15, —dE mag = —( E mag (t = tpeak + 10) — E mag (t = 0)) ~ 0.45, 
and dE int = E int (t = t peak + 10) — E int (t = 0) ~ 0.6. Note that the variation of the 
total energy remains within ~ 15 %. The field line rooted at ( x,y) = (1.8,0) also gives a 
similar result. The small variation in the total energy reflects the fact that the compressive 
and expanding motions by the surrounding plasma are localized in time and space (see 
Figure [9]). Therefore, the compressive motion of the reconnected flux tube is found to be 
not significant with respect to the total energy variation. 
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4. Pseudo-2D MHD Model of a Post-Flare Loop 
4.1. Physical Processes Considered 


Through the analysis of the 2D MHD simulation, we found that the field-aligned 
plasma motions (particularly evaporation flows and slow shocks) and heat conduction 
seem to mainly determine the dynamics in the post-flare loops. Fast shocks are important 
for converting the kinetic energy of the reconnection outflow to the heat and for locally 
changing the cross-sectional area, but they do not largely change the total energy of the 
field lines (Figure HU1). On the basis of the results, we aim to model the multi-dimensional 
reconnection and flare processes in a simplified MHD scheme. 


What is important for conduction-heating-type flare modeling would be 1. to model 
a reconnected field line, 2. to model the reconnection inflow and outflow, 3. to include 
MHD waves (since MHD waves can carry a large amount of the released energy from the 


reconnection sites flKigure et ah 


2010 )). 4. termination of the reconnection outflow and 


energy conversion of the kinetic energy of the outflow into the thermal energy, and 5. to 
include the heat conduction which is essential to determine the temperature of the flare 
loops. 


Considering that the plasma motions are frozen-in a magnetic field, a ID MHD model 
will be the simplest form among the possible MHD models. We developed a model based 
on the ID MHD equations, where all the variables depend on one space dimension and all 
the three components of the magnetic and velocity fields are considered. We regard our 
model as a pseudo-2D MHD model. Note that the meaning of ”pseudo-2D” of our model is 


different from that of 


Hori ct al. 


(119971) hydrodynamic model which consists of isolated and 


fixed semi-circular loops with different lengths and constant cross-section. 
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4.2. Assumptions and Basic Equations 

We take a cartesian coordinate system in which all the variables are functions of x. 
The reconnection outflow is in the y-direction. 0 is the ignorablc coordinate (uniform in the 
^-direction). 

We mimic a reconnected magnetic field line by assuming a magnetic field with a 
sharp bend (see Figure fTTl) . The magnetic field line shrinks with time and drives the flow 
perpendicular to the field which represents the reconnection outflow. 

A schematic picture of our model is as follows. Figure fl2l displays a comparison of our 
pseudo-2D MHD model with the 2D MHD model. h rx is the height of the reconnection 
point, and x max is the location of the foot-point of the field line. They are linked by the 
relation of x max = h rx tan 6. The reconnection angle 9 and the initial plasma beta /3 are 
treated as free parameters. If a guide-field (the 0 -component of the magnetic field) is 
included, the guide-field angle cj) = tan -1 (. B z /B x ) will be an additional parameter. 

The shrinking motion of the reconnected field line will stop when it collides with the 
magnetic loops piling up below. To model this process, the final configuration of the field 
line is given in the model and the outflow is decelerated when the field line approaches the 
final state. The termination process is modeled by adding a damping term to the equation 
of motion perpendicular to the moving field line. We let the damping term work only when 
the reconnected field line comes close to the final configuration. 

According to Figure flOl the total energy of a field line is conserved within several 10 %. 
On the basis of this result, we hypothesize that the total energy in a reconnected flux tube 
is conserved. The cross-sectional area in our model is assumed to be uniform and constant 


in time. 


The basic equations are as follows: 


de d 
dt dx 


P 


\B\ 


-)v x 


8n 4n 

9{pv, 


9p d x 
9B V 9 , 

~dt - ^ yBX - l ’ XBy) 

9Bz 9 

ry, rs \^Z^X ^X^z) 

at ox 

-Bjv • B) + k 0 T 5 / 2 B * 9T 


xj 9 

dt + dx 

9(pv y ) 

dt 

d(pv z ) 
dt 

e = 


\B \ 2 dx 
2 |-B | 2 5 2 ' 


s 


+ 


A 






V 


7 


B X By 

Ait 

B X B Z 

An 

1 , ,2 I - 5 ! 2 

T + 2 pM + ~sTT’ 

kB „ 
p = —pr, 
m 



(19) 


(20) 


(21) 


(22) 

u b (Ad(x,t))pv Xx , 

(23) 

u b (Ad(x,t))pv± y , 

(24) 


(25) 


(26) 


(27) 


where v = (v x ,v y ,v z ), B = (B x ,B y ,B z ), v± x = (v - U||) x , v ±y = (v - v\\) y , and 
V|| = B ( v ■ B)/\B\ 2 . All the physical quantities are only dependent of x and t. Note that 
we include a damping term in the momentum equations that slows down only the x and 
y-components of the velocity perpendicular to the magnetic field. The detailed description 
of the damping term is given below. Note that the total energy is conserved along a field 
line, and that the kinetic energy reduced by the damping term is converted into the thermal 
energy. The normalization units are the same as those in the 2D MHD model (See Tabled]). 
Note that this model can treat slow-mode, fast-mode and Alfven waves. 


Considering the symmetry, we only solve the domain within 0 < x < Xm ax . At the 
boundary of x = 0 the mirror symmetric boundary conditions are imposed. The boundary 
of x — x max is free. 


The initial conditions are as follows. The free parameters that determine the initial 
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magnetic field are the plasma beta 0 = 87rp/|i?| 2 , and the angles 6 and 0 (see Figure IT2]) . 
The guide held effect is detailed in Appendix [B] The initial magnetic held is given by 


B 0 

/8vrp 0 

V P ’ 

(28) 

B x (x) 

= — B 0 sin 6 cos 0, 

(29) 

By(x ) 

= B 0 cos 6 COS 0, 

(30) 

B z (x) 

= Bq sin 6 sin 0. 

(31) 


The gas pressure is assumed to be uniform (p 0 ). The domain is divided into the two layers, 
namely chromosphere and corona: 


P ( X ) Pcor T ( Pchr Poor) ^ 

where Xtr is the location so that yu n e( x TR,t = 0) = h TR . 
that describes the configuration of the held line. 


1 + tanh 


x - x T r 


(32) 


WTR 

Where yn ne (x,t ) is a function 


The initial and final magnetic held configurations are described as follows. See also 
Figure [131 The reconnection point (the location of the sharp bend) is assumed to be at 
(x, y) = (0, h rx ). The initial condition of a magnetic held is written as 


!/line(^d 0) 


tan 9 


lx 


X r . 




(33) 


for x > 0. The final state of the magnetic held is approximated by a quadratic function of 


Z/Zme,/m(*^) 


X r , 


2 tan 6 


x 


X n 


- 1 


(34) 


which has the same slope of the tangent line with the equation (|33|) at x = x r 


We virtually consider the travel distance of the reconnected held line in the ^-direction. 
When the held line approaches the final state, the damping term is applied only to the x- 
and y-components of the velocity perpendicular to the held line. We define the distance in 
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the ^-direction between the held line at (x, t) and the held line in the hnal state as 

Ad(:£, t) yiine(,%jt') (35) 


The damping term only works when the held line approaches to the hnal state: 


u d (Ad(x,t)) 


1 


td 


amp 


1 

2 


1 — tanh 


/ A d(x, t ) \ 

V w d )\ ’ 


(36) 


where tdamp = w d /vA , y , Va, v = B y /\J\^p cor is the outhow speed in the y-direction, and 
w d is a free parameter that denotes a typical braking distance. To prevent the held line 
from shrinking further even after the travel time t trave i = (h rx — yn n ejin{x = 0))/Va i3/ , we 
increase the damping coefficient v d by a factor of 100 after the time t = 1.2 ttravei- Note 
that the term arising from the damping terms is not included in the energy equation under 
the assumption that the total energy along a held line is conserved. The kinetic energy 
decreased by the damping term is converted only into the thermal energy. 


The reconnection angle 6 and plasma beta /3 are important parameters to determine 
the total released magnetic energy and energy conversion rate (say, reconnection rate). 
To cho ose a phys ically acceptable parameter set, we utilize the analytical approach by 


Fallc ct al. 


(119981) . A detailed description is presented in Appendix iBl 


The height of the reconnection point h rx is assumed to be 20 = y m ax, 2 Di and the 
reconnection angle 9 is 7 t/ 12. The domain size is therefore x max = h rx tan 9 = 20 tan (7 t/ 12). 
The width of the transition region and the typical damping distance are respectively 
wtr = u>tr, 2 D = 0.2 and w d = 2. 


The numerical scheme of the pseudo-2D MHD model is based on the 


Vogler et al. 


( 20051 ): the fourth-order space-centered difference for spatial derivative and an explicit 


four-step Runge-Kutta time integration. We explicitly solve the heat conduction term. 
Using the current computational resources, it is not difficult to explicitly solve the heat 
conduction in our ID calculations. The domain is resolved by 640 grids. 
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5. Numerical Results of Pseudo-2D MHD Model 
5.1. Dynamics and Energetics 

Figure fT4l demonstrates the time evolution of the magnetic held structure. The held 
line retracts and sweeps up the plasma like a slingshot. The reconnection outhow is 
decelerated when the held line approaches the assumed final configuration. In the following, 
we performed the same analysis as for the 2D MHD model. 

One-dimensional plots shown in Figure fl5l demonstrate the formation of a pair of two 
slow shocks attached to the reconnection outhow, that is, the Petschek-type slow shocks. 
As well as the 2D MHD model (Figure [3]) , the Petschek-type shocks are isothermal shocks 
due to the heat conduction. 

Figure [IB1 compares the field-aligned motions in the pseudo-2D MHD model and 2D 
MHD model. As well as in the 2D MHD model, the heat generated at the Petschek 
slow shocks is transferred to the chromosphere by the heat conduction, forming the 
chromospheric evaporation. The humps and high-density regions at the top are found to be 
formed in the same way as in the 2D MHD model: the downflow-evaporation collision and 
shock-shock interaction, respectively. The damping of the slow shocks are also observed. 

The energy evolution is examined in the pseudo-2D MHD model. When we integrate 
the energies, we drop off the term 1/B in the integrand because in the pseudo-2D MHD 
model the variation of the cross-sectional area is not considered. Figure [171 is the same as 
Figure fTTTl but for the pseudo-2D MHD model. It is shown that — dE mag ~ 0.6, dE int ~ 0.55, 
and Ekin^max 0.25, which is similar to the results of the 2D MHD model. 
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5.2. Dependence on Parameters 


Previously a formula that determines the post-flare loop temperature was derived 
under the assu mption that the energy input to a loop balances with the conduction cooling 


rate ([Fisher & Hawley 


1990). The formula is given by 

2/7 


QL 2 


loop 


V 


(37) 


2k 0 J 

where Q is the volumetric heating rate and Li oop is the half length of the magnetic field 
line of a post-flare loop. The heating rate by magnetic reconnection is determined by the 
Poynting flux: Q = B 2 /( 47 t ) x Va/Li oop . Using this, the temperature can be written as 


t « r 3/7 fcXV /7 , 


by assuming that Li oop ~ h rx . This scaling law was derived by 


(38) 


Yokovama fc Shibatal (1998 


200lh . We checked whether the scaling law based on a magnetic reconnection model holds 
in our pseudo-2D MHD model. 


Figure US shows the numerically-obtained f3-T, h rx -T and kq-T relations. Tj et denotes 
the temperature in the reconnection outflow. T* denotes the maximum temperature after 
the pair of the two slow shocks generated by the chromospheric evaporation flows collides 
at the apex. As shown in Figure [IS it is found that the temperature in the loop obeys the 
scaling law well. 


6. Summary and Discussion 

In this paper we investigated the flow structure, shock formation and thermal evolution 
in the post-flare loops using MHD simulations. On the basis of the 2D simulation result, we 
have developed a new post-flare loop model (the pseudo-2D MHD model). We compare the 
flow structure, shock formation and energetics of a specific held line in the 2D calculation 
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with those in the pseudo-2D MHD model, and then we find that they give similar results. 
Here we summarize the results and compare our results with previous studies. 


Performing a 2D MHD simulation, we found new shock structures. The termination 
shock consists of two oblique fast-mode shocks and sometimes a horizontal fast-mode 
shock (Figure S]). A hump is formed as a result of the collision of the downflow and the 
chromospheric evaporation flow. After the collision, the fronts of the evaporation flow and 
downflow become slow shocks, and then the hump appears as a dense region behind the two 
slow shocks (Figure [5]) . The upward slow shock finally interact with the slow shock coming 
from the other side at the top, forming the high-density region. Note that the high-density 


region is separated from the blob in 
(Figure [6]). 


Yoko vama fe Shi batal 020011 ) and is formed below it 


We found that the strength of the termination shock in the 2D MHD model shows 
a quasi-periodic oscillation, which is a multi-dimensional feature. In addition, the shock 
reflection and Mach reflection are sometimes observed in a concave magnetic structure, 
which could be important for understanding the heating in the loop-top. These complicated 
structures at the top will be detailed in our future papers. 


We observe no prominent shocks nor waves propagating back and forth from end 
to end. In addition, no prominent standing slow-mode waves are observed within the 
calculated time range. By performing 2D MHD simulations without the heat conduction, 
we confirmed that without the heat conduction the slow shocks formed at the fronts 
of the downflows from the top propagate back and forth from end to end. Thus the 
propagation of the slow shocks in the post-flare loops is found to be significantly affected 
by the heat conduction and evaporation flows. The slow shocks are damped by the heat 
conduction. The propagation speed of the slow shocks is reduced by the evaporation flow 
(Doppler effect), which makes the shock propagation time longer. Therefore it is essential 
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to consider the flows resulting from reconnection (particularly downflows from the top) for 
understanding the behavior of magneto-acoustic waves in the post-flare loops. 


In fully 3D situations, reconnection can be intermittent in space and time, which 
could affect the shock structures found in this study. 3D component reconnection, 
where reconnecting magnetic field lines are not perfectly anti-parallel, can result in the 
reconnection outflow jet with a speed insufficient for the formation of the fast shock above 
the loop top. However, in the case that the reconnection outflow speed exceeds the fast 
mode phase speed in the outflow region, we expect the formation of the termination shock 
found as in our 2D MHD model. We also expect that the slow shocks presented in the 2D 
MHD model will be formed in 3D, since the field-aligned plasma flows which are essential 
to form slow shocks are well treated in our MHD flare modeling scheme. 


From the 2D simulation, we found that the field-aligned plasma motions (particularly 
evaporation flows and slow shocks) and heat conduction mainly determine the dynamics in 
the post-flare loops. Considering this, we construct the pseudo-2D MHD model which is 
basically a ID MHD model. The pseudo-2D MHD model is compared with the 2D MHD 
model, and we found that the dynamics (particularly flow structure and shock formation) 
and the energetics are similar between the two models. The scaling law for the temperature 
based on a reconnection model is also examined, and it is found that the scaling law holds in 
the pseudo-2D MHD model. These facts indicate that our pseudo-2D MHD model captures 
important features of a reconnection model of solar flares. 


ID hydrodynamic models, like 


Mariska et ah 


( 19891 ) and 


Hori ct ah 


(1997 


1998!), have 


been used for the post-flare loop modeling. The models are useful to study the thermal 
evolution and flows in the flare loops, but the energy input (in many cases the heat input) 
must be done by hand. Our model includes many features of the multi-dimensional MHD 
processes related to magnetic reconnection, like the heating by the Petschek slow shocks and 
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conversion of the kinetic energy of the reconnection outflow to the heat. Another important 
point is that our model can treat MHD waves and shocks generated in the post-flare loops, 
which could be important to understand the flow structure in the loops. 


We note that all the previous ID hydrodynamic models lack the strong downflows from 


the top (for example, see Figure 5 in 


Hprietah 


(119971) ). We showed that the downflows play 


important roles in forming the humps and slow shocks (Figures [5] and [6]). Our pseudo-2D 
MHD model naturally produces the downflows from the top, allowing us to study the 
dynamic evolution of the thermal structure. 


A theoretical model in which held lines shorten a: 
on the thin flux tube approximation was proposed by 


'ter localized 3D reconnection based 


Longcope et ah 


( 20091 ). The model 


assumes that the plasma has always low-/?. Our 2D MHD simulation demonstrates that, 
because of the shock heating and compression, the plasma /? becomes larger than unity 
not only in the reconnection outflow but also in a large part of the post-hare loops (see 
Figure [7]). In such regions, we need to consider the effects of the gas pressure to understand 
the how and shock structures. Our pseudo-2D MHD model can treat the high-/? plasma. 


Lon gco p e et ah 


(120091) model could be valid in the situation in which the shock heating 


does not break the low-/? assumption (e.g. in the situation in which the guide-held is much 
stronger than the reconnection held). Regarding this issue, see also Appendix iBl 


The high-/? condition could lead to disordering of the magnetic held when turbulence 
is important so that the geometrical assumptions made for our pseudo-2D MHD model are 
not met. We found no evidence that turbulence is important in the current sheet in our 
2D MHD model. Our model will be useful to model the reconnected held lines in such a 
situation in which turbulence is not important. 


As a result of the approximations made for the pseudo-2D MHD model, some 
multi-dimensional processes such as termination shock formation and turbulence cannot 
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be modeled. However, the pseudo-2D MHD model are able to treat plasma flows and 
waves/shocks (slow-mode waves/shocks and Alfven waves) along the magnetic field, which 
conld be useful to understand the plasma motions in the post-flare loops. 

Our pseudo-2D MHD model requires much smaller computational cost than other 
multi-dimensional MHD models. This model will allow us to study the evolution of the 
flare loops in a wide parameter space without expensive computational cost. Also, it will 
include detailed physics like the non-equilibrium ionization effect (e.g. 

. These will be our future work. 


The authors are grateful to Dr. Seiji Zenitani for fruitful comments on shocks found in 
our simulations. ST acknowledges support by the Research Fellowship of the Japan Society 
for the Promotion of Science (JSPS). This work was supported by a Grant-in-Aid from the 
Ministry of Education, Culture, Sports, Science and Technology of Japan (No. 25287039). 

Facilities: 


be much easier to 


Irnada et ah, 


2011) 


A. Numerical Method of Heat Conduction in 2D MHD Model 


We modify the tim e st ep splitting method to calculate the heat conduction flux by 


Yokovama fc Shibata ( 200li ). Using the MHD energy flux F m hd and heat conduction energy 


flux F c , we can write the energy equation as follows: 

V • F c = 0, 


dE „ „ 

o, + V • 

at 


(Al) 


where E is the total energy. The discretized form is 


— (E n+1 - E n ) + (V • F mhd ) n+1 / 2 + (V • F c ) n+1/2 = 


0 . 


(A2) 
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First, we calculate the MHD part 

-Uf* 

a v 


E r 


= -(V-F, 


mhd . 


i 71+1/2 


(A3) 


where the superscript * denotes the results of the MHD step. Then we calculate the heat 
conduction step 

1 


At 


(. E n+1 - E*) = -(V ■ F c ) n+1/2 . 


(A4) 


The procedure mentioned above is the same as 


Yokova m a fe Shibata ( 2001 ). The 


differen ce appears in t 


Yokoyama & Shibata (2001) are 


ie_expression of the heat flux formula. The heat flux formulae in 


73 * 73 * 

- Ko (r*) 5 / 2 —(— ■ vt), 


ST 

_A* _ _ A* 

xx dx xy d y ’ 


F, 


c,y 


-A 


yy 


dT 

dy 


A 


yx 


dT 
' dy ’ 
dT 
dx ’ 


r>* r>* 

A* _ (r r *\5/2 JD a JD b 

^ab K 0D ) (B *) 2 ’ 


(AS) 

(A6) 

(AT) 
(AS) 


where the subscripts a and b denote x and y. The time and space discretization are, for 
example, 


dT dT 

A* t U± V^+l A* f U± 
~ AxA ^ ] ~ A 


xy\ 


dy 


We modify this to 


F r 


dT dT 

_ a * _ a * /y__'\ n +i 

xxK dx’ xyK dy’ ' 


(A9) 


(A10) 


The difference of two formulae appears in the operator matrix of the implicit scheme 
used. The number of the non-zero components in one row is 5 in Equation (lA9jh and 9 
in Equation (IA10I) . Geometrically, in Equation (lAlOj) we use all the neighboring 8 grids 


around each point at (n + 1) step to calculate the heat flux (Figure [19]). This method gives 
more accurate results than the previous method, particularly in the grid points where a 
magnetic field is largely bend and oblique to coordinate. 
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B. Analytical Approach to the Riemann Problem: Possible Solutions 

We shall describe the exact solutions of the symmetric MHD Riemann problems that 
can appear in our pseudo-2D MHD model. In a symmetric MHD Riemann problem, the 
initial states in the left and right hand side are assumed to be as follows. 


PL 

— Pr, 

(Bl) 

PL 

= PR, 

(B2) 

B X ;L 

B x .r, 

(B3) 

By-L 

= ~ By-R , 

(B4) 

B z -l 

= b z . r , 

(B5) 

Vfi 

= v L = 0, 

(B6) 


where the subscripts L and R represent the left hand side and the right hand side, 
respectively. 

Let us consider the situation without the guide held component ( B Z . L = B z . R ~ 0)- 
Since the magnetic held vectors in the initial states are in the x-y plane, the rotational 
discontinuities will not appear in the solutions. Due to the initial discontinuity, one contact 
discontinuity, two slow shocks, and two fast rarefaction waves are generated. 

From the jump condition across the contact discontinuity, the tangential magnetic 
held, By, should be equal to zero in the region between two slow shocks. Therefore, it is 
necessary to have either the switch-off slow shock (SS) or the switch-off fast rarefaction 
wave (FRW). In the former case (see the left panels of Figure 120]) . the exact solution consists 
of the switch-off SS and FRW. On the other hand, in the latter case (the middle panels of 
Figure I2U1). the exact solution consists of the pure hydro shock (HS) and switch-off FRW. 

There are two parameters in the symmetric MHD Riemann problems: plasma beta 
/3 and the reconnection angle 6. A sophisticated procedure to derive the exact solutions 
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of more general MHD Riemann problems was developed by 


Falle et al 


(1998). We utilize 


the procedure to derive the exact solutions of our problems. The result is summarized in 
the phase diagram of Figure EH where f3 and 9 are the key parameters. Here the heat 
conduction is neglected. The 9-/3 space is divided into two regions: switch-off slow shock 
regime and pure hydro shock regime. The parameter set of the example introduced in this 
paper is within the switch-off slow shock region. If we include the heat conduction effect, 
the pure hydro shock region will become slightly wider. 


Figure EH displays the normalized reconnection rate in the 9-/3 diagram, where the 
normalized reconnection rate is defined by v y B x /(V aB x0 ). Where v y is the outflow speed, 
Va = y/B% + B/; /y/p, and B x = B x0 . The maximum reconnection rate observed in solar 
flares and simulations is of the order of 0.1. To take the reconnection rate similar to 0.1 in 
our model, we need to choose a reasonable parameter set from this diagram. 


If a guide held exists, the rotational discontinuities (RD) will appear in addition to 
SS and FRW. In this case it is necessary to have either the switch-off SS, switch-off FRW, 
or RD because of the boundary condition at the contact discontinuity ( B y = 0). Among 
them, only the solution with RD (the right panels of Figure [20l) meets the requirement 
that the continuity condi ti on of V x a t the contact discontinuity should be satisfied 


(Pets chek & Thorne 1967). 


Longcopa ( 20141 ) claimed that pure hydro shocks can be formed 


in a reconnected flux tube even when a guide held exists. However, considering the analysis 
here, pure hydro shocks will not appear when a guide held exists. 
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Tabic 1: Normalization units 


Quantity 

Unit 

Value 

Length 

L 0 

3,000 km 

Velocity 

c s0 = [7 (k B /m)To] 1/2 

170 km s -1 

Time 

Lq/ C s o 

18 s 

Temperature 

To = T cr 

2 x 10 6 K 

Density 

Po Per 

( 10 9 cm ' 3 )xm 

Pressure 

'y(k B /m)poT 0 

0.47 dyn enr 2 

Magnetic field strength 

[8TT'y(k B /m)p 0 T 0 ) 1/2 

3.4 G 


® z 


Y 


Hot tenuous layer 
(Corona) 

Cool dense layer 
(Chromosphere) 

Fig. 1.— Schematic diagram of the initial and boundary conditions of the 2D MHD simu¬ 
lation. The solid arrows denote the magnetic field. The circled region indicate the region 
where the localized resistivity is applied. The grey region at the bottom is the chromosphere. 
The calculated domain is the region where x > 0. 



Free boundary 
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Temperature 

t=12.0 t=13.2 t=15.0 t=20.0 



Weak disturbance 
due to the initial condition 

14 


Log Density 


6 

Chrom. 

evaporation 2 

o 


If W m m 

1 j f j \\\ \ / ///J \ \\\V/// 

k\\ 

1 

y 


-4 -2 0 



Fig. 2.— Time evolution of the temperature (Top) and density (Bottom) of the 2D MHD 
model. The solid lines denote the magnetic held. 
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One-dimensional plots of the reconnection outflow 

(2D MHD Model) 

Time = 13.0, Y= 10.0 



Fig. 3.— One-dimensional plots parallel to the x-axis across y — 10 at t — 13 (across the 
reconnection outflow region). The density, temperature, pressure, J z , v x and v y are displayed. 
The vertical dashed-lines indicate the slow shocks attached to the reconnection outflow. 
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div V/C s 
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- 2-1012 

X 
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Oblique 

fast-mode shocks 



13 

12 


fast-mode shock (Mach disk) 


Oblique 

fast-mode shocks 

1 2 


Fig. 4.— Gas pressure (Left) and V • v/C s (Right) maps of the loop top. Top: Snapshot 
when the termination shock consists of two oblique fast-mode shocks. Bottom: Snapshot 
when the termination shock consists of two oblique shocks and a horizontal shock. Note that 
most of the regions where V ■ v/C s takes large negative values are slow or fast shocks. 
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Formation of hump 


Logio Density div V/Cs 


- 0.5 0.0 0.5 1.0 
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Slow shock > 
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o 

X 


1 
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Front of downflow 

Front of chrom. 
evaporation 


Hump 


Slow shock 


Fig. 5.— Formation of the high-density regions called humps. Left: the density in loga¬ 
rithmic scale. Right: V • v/C s . The solid lines denote magnetic held lines. The velocity 
vectors projected to the x-y plane are also displayed in the density maps (only the vectors 
with yjv* + Vy > 0.3 are shown). Note that most of the regions where V ■ v/C s takes large 
negative values are slow or fast shocks. 
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Formation of high-density loop-top 
Logio Density div V/Cs 


- 0.5 0.0 0.5 1.0 




Oblique fast shock- 
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Oblique fast shock 


Shock reflection 


High-density region 


Fig. 6.— Formation of the high-density loop-top. Left: the density in logarithmic scale. 
Right: V • v/C s . The solid lines denote magnetic held lines. Only one of the pair of the slow 
shocks heading to the apex is indicated by the arrows. Note that most of the regions where 
V ■ v/C s takes large negative values are slow or fast shocks. 
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Log 10 Density 


- 0.5 0.0 0.5 1.0 



Temperature 
0.5 2.0 3.5 5.0 



- 4-2 0 2 

Entropy 

-2 -1 0 


- 4-2 0 2 

Log 10 P 
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Fig. 7.— Snapshots of the density, temperature, entropy (log (p/p 7 )), and plasma (3 at 
t = 20. The arrows indicate slow shocks in the post-flare loops (See also Figures [5] and |6|). 
The solid lines in the first three panels denote the magnetic held. The contour in the plasma 
f3 map denotes the level where (3 = 1 (log 10 (3 = 0). 
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Logio Density 

Time = 10.0 Time = 12.0 Time = 14.0 Time = 15.0 Time = 20.0 



01 2301 2301 2301 2301 23 

X X X X X 


Fig. 8.— Evolution of a specific reconnected field line with the density map. The tracked 
field line is indicated by the thick black line. The foot-point of the field line is located at 
(x,y) = (1.2, 0). The circle is the origin of the distance. 
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2D MHD Model 


Log, 0 Density 

- 0.5 0.4 1.3 


Temperature 

0 2 4 



High-density region 
(shock-shock interaction) 


Hump 

(downflow-evaporation collision) 

Reconnection 
outflow region 

Disturbance — 
due to 1C 


parailel 



Temperature increase 
due to the thermalization of 
the kinetic energy 

Heat conduction front 
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Chromospheric 

evaporation 


-1 0 


-1 0 1 



Slow shock 


Slow shocks formed 

after downflow-evaporation collision 

Chromospheric 

evaporation 

Petschek shock 

Disturbance 
due to 1C 


Fig. 9.— Time-distance diagrams of the density (Top Left), temperature (Top Right), uy 
(Bottom left), and V • v/C s (Bottom Right). The physical quantities are measured along 
the held line which is originated from (x,y) = (1.2,0). The distance is measured along the 
held line, and its origin is the apex (x — 0). The data of the domain of x < 0 is also shown 
just for visual inspection. The sign of uy is defined as positive and negative when a plasma 
how travels to the apex and foot-points, respectively. This held line does not pass through 
the termination shock. Note that the regions where V • v/C s takes large negative values are 
slow or fast shocks. 
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2D MHD Model 

Xfoot =1.2 



Time 


Fig. 10.— Time evolution of the total (thick solid), magnetic (dashed), internal (dash-dot), 
and kinetic (thin solid) energies along the specified held lines in the 2D MHD model. The 
vertical dashed lines denote the times when each held line reconnects. 
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2D picture 1D approximation 
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Fig. 11.— Schematic diagram of the reconnected magnetic held in 2D (Left) and the mag¬ 
netic held mimicking the reconnected held in the lD-approximation model (pseudo-2D MHD 
model, Right). The bottom panels show the difference in the reconnected magnetic held con¬ 
figurations between two models. 
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2D picture Pseudo-2D MHD model 




Fig. 12.— Schematic diagram of the set up of our pseudo-2D flare model (Right) based on 
a 2D picture (Left). h rx and hrn are the height of the reconnection point and height of the 
transition region, respectively, tan -1 (x max /h rx ) = 9 is the reconnection angle. The grey 
shaded regions indicate the dense cool material representing the chromosphere. Note that 
in the pseudo-2D MHD model, all the physical quantities are function of x. 
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Y 



Fig. 13.— The magnetic configurations of a specific field line in the initial state (yu n e(x,t 
0 )), at t — t ( yu n e{x,t = t )), and in the final state ( ynnejin(x )), respectively. A d(x,t) 
yiine(x,t ) - yiine,fin(x) is also described. 
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X 


Fig. 14.— Downward moving magnetic field line in the pseudo-2D MHD model. The dashed 
line denotes the assumed final configuration yu ne ji n ■ 
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One-dimensional plots of the reconnection outflow 
(Pseudo-2D MHD Model) 
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Fig. 15.— One-dimensional plots parallel to the x-axis at t — 0.7 (across the reconnection 
outflow region). The density, temperature, pressure, J Z) v x and v y are displayed. 
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Fig. 16.— Comparison between the pseudo-2D MHD model (Left) and the 2D MHD model 
(Right). Time-distance diagrams of the density and V • v/C s are shown. The physical 
quantities of the 2D MHD model are measured along the held line which is originated from 
(x,y) = (1.2,0). The distance is measured along the held line, and its origin is the apex 
(x — 0). The data of the domain of x < 0 is also shown just for visual inspection. Note that 
most of the regions where V ■ v/C s takes large negative values are slow or fast shocks. 
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Pseudo-2D MHD Model 
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Fig. 17.— Time evolution of the total (thick solid), magnetic (dashed), internal (dash-dot), 
and kinetic (thin solid) energies along a field line in the pseudo-2D MHD model. 
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Pseudo-2D MHD Model 



10 

Log 10 h rx 



Fig. 18.— Temperature as a function of the initial (3 (Top), the height of the reconnection 
point h rx (Middle), and the heat conductivity k 0 (Bottom), respectively. The data points 
are taken from the pseudo-2D MHD simulations. T ]et denotes the temperature in the re¬ 
connection outflow. T* denotes the maximum temperature after the pair of the two slow 
shocks generated by the chromospheric evaporation flows collides at the apex. The dashed 
lines show the analytical scaling laws. 
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Yokoyama & Shibata 2001 
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Our method 



• Points with the values at t = t* 

• Points with the values at t = t n+1 


Fig. 19.— The grid points used to calculate the heat flux for the grid (i,j) in 


Yokoy ama fe Sh i batal (120011 ) (Left) and our method (Right). The red and blue grids have 
the physical quantities at the time t = t n+1 and those at t — t*, respectively. 
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Without a guide field 
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Fig. 20.— Three possible solutions of the Riemann problem. Left column: Case without a 
guide field, where B y is switched off by the slow shock (SS). Middle column: Case without a 
guide field, where B y is switched off by the fast mode rarefaction wave (FRW). Right column: 
Case with a guide field, where the magnetic field is rotated by the rotational discontinuity 
(RD) and the magnetic energy is released at the slow shock (SS). First row: x — t diagrams 
that show the Riemann fans. Second row: \B Z \ — \B y \ diagrams that show the transitions 
through the regions in the Riemann fans (the paths are indicated by the dashed lines in 
the first row). The third and fourth rows: x — y diagrams that show the magnetic field 
configurations on the plane. Only for the case with a guide field, x — z diagram is also 
depicted. The grey regions indicate the high-entropy regions due to the shock dissipation. 
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Fig. 21.— A phase diagram of the possible solutions of the symmetric Riemann problem in 
the situation without a guide-held. 6 is the reconnection angle and (3 is the plasma beta. 
”Switch-off slow shock” regime and ’’pure hydro shock” regime corresponds to the left and 
middle panels in Figure [20] respectively. 
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Reconnection Rate 




reconnection rate for (3 = 0.1 as a f un ction of 6 . 
















